Quantum Ratchet Accelerator without a Bichromatic Lattice Potential 
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In a quantum ratchet accelerator system, a linearly increasing directed current can be dynamically 
generated without using a biased field. Generic quantum ratchet acceleration with full classical chaos 
[Gong and Brumer, Phys. Rev. Lett. 97, 240602 (2006)] constitutes a new element of quantum 
chaos and an interesting violation of a sum rule of classical ratchet transport. Here we propose 
a simple quantum ratchet accelerator model that can also generate linearly increasing quantum 
current with full classical chaos. This new model does not require a bichromatic lattice potential. 
It is based on a variant of an on-resonance kicked-rotor system, periodically kicked by two optical 
lattice potentials of the same lattice constant, but with unequal amplitudes and a fixed phase shift 
between them. The dependence of the ratchet current acceleration rate on the system parameters 
is studied in detail. The cold-atom version of our new quantum ratchet accelerator model should 
be realizable by introducing slight modifications to current cold-atom experiments. 

PACS numbers: 05.45.Mt, 05.45.-a, 05.60.Gg, 32.80.Qk 



I. INTRODUCTION 

A ratchet accelerator (RA) [l| can generate, with- 
out using a biased field, directed transport in both the 
momentum and coordinate space. Specifically, certain 
spatio-temporal symmetries in the Hamiltonian dynam- 
ics are broken and as a result a linearly increasing di- 
rected current can be dynamically generated. Such a 
property of RA systems is of considerable interest for un- 
derstanding (i) general properties of quantum and classi- 
cal ratchet effects in Hamiltonian systems [11, d, H, 0, H, H, 
S H, P| , (ii) quantum-classical correspondence in trans- 
port phenomena and (iii) a number of interesting 
topics in quantum chaos p^ . 

Ongoing cold-atom studies of the well-known quantum 
kicked-rotor (QKR) model [iq| have motivated several 
RA studies using QKR variants. In particular, Ref. 
showed that accelerating quantum ratchet current can be 
realized by considering a QKR variant, with the kicking 
period on the main resonance with the recoil frequency 
of the cold atoms. Reference showed that a quantum 
RA can also be realized with QKR variants on high-order 
quantum resonances. In both cases, spatio-temporal 
symmetries in the dynamics are broken by using a bichro- 
matic optical lattice (specifically, an optical super-lattice 
obtained by superimposing two stand-waves with peri- 
ods A/2 and A/4). However, though already achieved in 
some static cases [ll|, experimentally realizing bichro- 
matic and pulsed optical lattices is somewhat demand- 
ing. Indeed, in two recent cold-atom on-resonance- QKR 
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experiments |12l . [l3j . directed quantum transport is ef- 
fectively demonstrated by use of a single-period optical 
lattice only, with the price being that special symmetry- 
breaking initial superposition states should be prepared. 

The ratchet transport in the above-mentioned on- 
resonance-QKR models 0, H, [H, [H, occurs only for 
isolated values of the effective Planck constant (to be 
defined below). By contrast, using variants of another 
paradigm of quantum chaos, namely, the kickcd-Harpcr 
model [l5|. Gong and Brumer [l^ proposed a quantum 
RA model that works for an arbitrary value of the ef- 
fective Planck constant. In this sense, the ratchet trans- 
port in this new model [T6l| is generic. Furthermore, this 
generic RA model works even when the underlying clas- 
sical dynamics is fully chaotic, a situation where classi- 
cal ratchet transport necessarily vanishes according to a 
classical "sum rule" 0, [E] . Hence the work in Ref. [T^ 
represents an interesting and generic quantum violation 
of a classical theorem. 

The detailed aspects of the above-mentioned quantum 
violation of the classical sum rule are yet to be explored. 
Along this direction, a cold-atom realization of a generic 
quantum RA model would be of great interest. Never- 
theless, such experiments were thought to be challenging 
because the model proposed in Ref. [lB| also employed 
a hashing bichromatic optical lattice and it was unclear 
how a kicked-Harper-like model can be realized in a cold- 
atom laboratory. 

Thanks to our recent finding [TtI . [isj that exposed a 
direct connection between QKR and a class of kicked- 
Harper-like models, here we are able to propose a quan- 
tum RA model that (i) contains all the important in- 
gredients as the model proposed in Ref. [1^1, (h) does 
not require a bichromatic lattice potential, and (iii) is re- 
alizable by slightly modifying existing cold-atom exper- 
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iments of QKR dynamics. Indeed, this new RA model 
only requires an on-resonancc variant of QKR, kicked 
by two optical lattice potentials of the same lattice con- 
stant, but with unequal kicking amplitudes and a fixed 
phase shift between them. In addition to offering a sim- 
pler quantum RA model that is of theoretical interest, it 
is hoped that our results below will motivate cold-atom 
experimental studies in the near future. 

This paper is organized as follows. In Sec. II we show 
how a wide class of twisted kicked Harper models can be 
realized by using an on-resonance "double-kicked" rotor 
model. General discussions in Sec. II directly lead to 
an atom-optics proposal for realizing the RA model pro- 
posed in Ref. tlfl]. In Sec. Ill we simplify the cold-atom 
RA realization in Sec. II, resulting in a new RA model 
that does not need a bichromatic optical lattice potential. 
We then present and discuss detailed numerical results of 
our new RA model, with an emphasis placed on the de- 
pendence of the current acceleration rate on the system 
parameters. In Sec. IV we briefly discuss one extension 
of this study. Section V concludes this work. 



II. COLD-ATOM REALIZATIONS OF A WIDE 
CLASS OF TWISTED KICKED HARPER 
MODELS 

Our starting point is the so-called double kicked-rotor 
model (DKRM) [H, HO, [Ml that has been experimen- 
tally realized. We use scaled and dimensionless variables 
throughout. The DKRM Hamiltonian is then given by 



DKRM 



^+KVKiq)Y.6it-nT) 

n 

LVLiq)Y,^it-^T-ij), 



(1) 



where q (g [0, 27r)) and p are conjugate coordinate and 
momentum operators, T is the period for both kicking 
sequences, rj is the time delay between the two kicking 
sequences, K and L characterize the amplitudes of the 
kicking fields, Vr-(9) and Vl('?) are periodic functions 
of q with the period 27r. The associated quantum map 
C/dkrm for a period from nT + 0~ to {n + 1)T + Q~ is 
given by 



C/dkrm 



(2) 



where h represents an effective and dimensionless Planck 
constant for the DKRM system (hence p = —iTid/dq). 

Theoretically, we shall first focus on an ideal situation 
where cold atoms are injected with exactly zero quasi- 
momentum (2^ . With that simplification we may con- 
sider only a Hilbert space satisfying the periodic bound- 
ary condition associated with q ^ q-\- 27r. The quantum 
resonance condition Tfi — 47r then leads to 



(3) 



reducing J/dkrm to C/dkrM' 



(4) 



where we have defined the rescaled momentum 

P = VP (5) 



and the rescaled kicking amplitudes 

k = i^K, 
L = TjL. 



(6) 
(7) 



Due to the above momentum rescaling, the effective 
Planck constant now becomes 



h — rjfi. 



(8) 



Equations (|6l7l8p show that the rescaled dimension- 
less system parameters K, L, and fi can be easily tuned 
by adjusting the time delay between the two kickin g se - 
quences. Based on a previous DKRM experiment p^ . 
we estimate that in experiments the kicking amplitudes 
K and L can vary in the range of 0.1 — 100, and the effec- 
tive Planck constant Ti can at least vary in the range of 
0.05 — 20. Our computational studies in the next section 
will be based on these two ranges. 

To gain insights into the quantum-resonance-reduced 
quantum map in Eq. ([4]), let us first rc-intcrpreted it as 
follows. Reading the four factors in Eq. Q from right 
to left, one sees that within each period T, in effect the 
system is first subject to one kick, followed by a free 
evolution of duration unity; then the system is kicked a 
second time, followed by a second free evolution of the 
same duration, but now with the free Hamiltonian given 
by 



fffrcc = -pV2. 



(9) 



Such an effective Hamiltonian with a negative kinetic en- 
ergy term was first considered in Ref. [1^. With this 
interpretation, one may define an "?7-classical" limit of 
this quantum map, i.e., the h = rjh ^ limit with fixed 
K and L. This terminology is inspired by the so-called 
"e-classical" limit in early studies of QKR models in the 
presence of gravity [2^. Let q'^ and be the counter- 
parts of q and p in this "ry-classical" limit, with their val- 
ues right before t — nT denoted by and p^^. Further 
defining 



P 



(10) 



one easily finds the classical map associated with the "77- 
classical" limit, 



dqn 



In+l 



(11) 

(12) 
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In terms of the canonical pair g'^ and V^, the classical 
Hamiltonian i/,^ that generates this "ry-classical" map is 
then given by 

n 

Consider now the simplest choice for the kicking 
potentials, i.e., VK{q) — VL{q) — cos{q). Such a 
choice under the restriction K — L was adopted by 
the original experiment ^19] an d p revious theoretical 
studies of off-resonance DKRM do, Substituting 
^K^q) = Vhiq) — cos(g) into Eq. ([13]), the resulting 
"7y-classical" Hamiltonian becomes precisely the classical 
kicked Harper model in terms of and q'^ [l7| . Return- 
ing to the old representation the obtained kicked 
Harper Hamiltonian becomes 

H^j^jj ^ L cosiP" + q'^) + k cos{q'')^5{t - n). (14) 

n 

Comparing the Hamiltonian i?xKH with the standard 
kicked-Harper Hamiltonian as a function of q"^ and p*^, 
i.e., 

= Zcos(p^) +Xcos(g^)^(5(< - n), (15) 

n 

we can regard i?xKH ^ "twisted" version of the stan- 
dard kicked-Harper model H^^. With this in mind, 
the quantum map in Eq. Q in the case of VK{q) — 
VL{q) — cos(g) can be regarded as a quantized version of 
the twisted kicked Harper model. 

We now apply this on-resonance DKRM strategy to 
realize a twisted version of the quantum RA model pro- 
posed in Ref. [lB|- This RA model involves the classical 
Hamiltonian 

■^BKH Icosip") 

+K [cos(g'= + 4>i)+ sin(2g'= + <l>2)]^S{t - n). 

n 

(16) 

If we now consider the following scenario: 
Scenario I: 

VL{q) = cos{q); (17) 
VK{q) = [cos(q + (^i) + sm{2q + (/-s)], (18) 

then a twisted version of ^^bkh (i-^-j in terms of T"^ and g"^ 
rather than p'^ and g*^) naturally emerges from Eq. (jlSp . 
Clearly then, at least for a twisted version, a cold-atom 
quantum version of the bichromatic generalized kicked 
Harper model -ffgKH realizable, provided that a kicking 
bichromatic lattice potential such as [cos(g-|-<^i)-f sin(2(7+ 
(j>2)] J2n^i^ ~ '^^^ ^® realized. In the next section a 
simpler realization of quantum RA is obtained. 

Before ending this section we make one important re- 
mark. In the standard kicked Harper model in 



Eq. (jl5p . the momentum variable is an abstract canon- 
ical variable. This becomes obvious if we consider the 
canonical equations of motion, yielding that the moving 
speed in the coordinate space is not proportional to the 
momentum. As such, it is unclear whether the momen- 
tum variable in the kicked Harper model can be directly 
related to the mechanical momentum of a moving parti- 
cle. Dana managed to connect this abstract momentum 
variable with the mechanical momentum of a charged 
particle kicked by a special sequence of magnetic fields 
|25| . Here, through the cold-atom realization of a wide 
class of twisted kicked Harper models, we are linking the 
momentum variable in the kicked Harper model with the 
mechanical momentum of cold atoms. Only through such 
connections can the expectation value of the momentum 
be interpreted as a current of moving particles. 

III. RATCHET ACCELERATOR WITHOUT A 
BICHROMATIC OPTICAL LATTICE 

Our discussions in the previous section make it clear 
that, in realizing a wide class of kicked-Harper-like mod- 
els with on-resonance DKRM, the following canonical 
transformation or twist is necessarily involved: 

(9^p^)-(9^n• (19) 

Due to this phase space twist, the resultant systems 
should assume different symmetry properties than those 
analyzed in terms of q'^ and p'^. Hence, the symmetry- 
breaking considerations in Ref. no longer apply to 
twisted kicked-Harper-like models. As a result, the use 
of a bichromatic optical lattice as in Ref. may not be 
the simplest approach for symmetry-breaking. It is this 
recognition that motivated us to seek a realization of a 
quantum RA without using a bichromatic lattice poten- 
tial. This attempt is also consistent with a recent study 
of ratchet transport (in coordinate space only) using an 
off-resonance DKRM involving two optical lattices of the 
same lattice constant [26j . 

Specifically, here we shall demonstrate that an on- 
resonance DKRM with the following alternative scenario. 

Scenario II: 
VKiq) = cosiq); (20) 
yL(g) =cos(g + 0) (21) 

can already give rise to a simple quantum RA model if 
K ^ L. In addition to the on-resonance condition, this 
scenario only needs to introduce two small modifications 
to a previous DKRM experiment [l^. First, the two 
kicking sequences of optical lattice potentials should have 
different amplitudes. Second, there should be a fixed 
phase shift <j) between these two optical lattice potentials. 

Using Eq. ([4]), scenario II described above gives the 
following quantum map, 

UnA = e'^e-'^''°'^''+'^^e-'^e-'f'^°''(i\ (22) 
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Using Eq. one then obtains the 77-classical Hamil- 

tonian of this quantum map, 

i/j^A ^ L cosiq" +p'' + (l}) + k cos(g'=) ^ S(t - n) . (23) 



Let \n) be the eigenstates of the momentum operator 
p, with an eigenvalue nh for a Hilbert space with the peri- 
odic boundary condition. In the foUowing we shaU focus 
on the RA dynamics for the initial state |0). The classi- 
cal analog of this initial state is a classical ensemble with 
p'^ = and a random uniform distribution of q'^. Such 
an initial state is a trivial state because it is symmetric 
upon time-reversal operations or space-reflection opera- 
tions. With this choice of the initial state, any induced 
current afterwards must be due to some broken spatio- 
temporal symmetries in the ensuing dynamics. It is also 
worth noting that the system described by Eq. ([221) is 
invariant under the transformations q 2Tr — q, p —p 
and 4> — > — (/>. As a result the current should undergo a 
sign change under </> —> —0, leading to the expectation 
that the generation of a ratchet current is forbidden for 
= or (/) = TT. For this reason we focus on other values 
of (h. 




FIG. 1: Time dependence of the ratchet current for both the 
quantum ratchet accelerator model of Eq. (|22|l (solid lines) 
and its r;-classical limit (dash-dotted lines) described by Eq. 
(|23p . For panel (a), system parameters are h — 1, K = 3, 
L = 1, and (fi = n/2 (upper two curves), 7r/3 (bottom two 
curves). For panel (b), system parameters are h — 1, K = 4, 
L = 2, and = n/2. The classical current in panel (b) 
remains indistinguishable from zero at all times because the 
system is in the full chaos regime (see Fig. 3d). 



A. Examples of Accelerating Ratchet Current 

Consider first a few computational examples depicted 
in Fig. 1. There, the time dependence of the quantum 
ratchet current, i.e., the expectation value of the scaled 
momentum, (p), is shown for some particular values of K, 
L, and (p. The initial state is |0) and the time propagation 
due to the quantum map C/ra is carried out by standard 
fast-Fourier-transformation techniques. The cases shown 
in Fig. 1 display spectacular linear acceleration of the 
ratchet current at a significant rate. In order to have a 
comparison with the underlying ry-classical limit, we have 
also calculated the ensemble-averaged classical momen- 
tum, {p'^). The classical calculations are based on the 
77-classical map given by Eq. (j23p . using an ensemble of 
10^ particles initially distributed along p'^ = randomly 
and uniformly. As is seen from Fig. la, the ry-classical 
currents can also increase linearly, with a slope smaller 
than their quantum counterparts. The result in Fig. lb 
is even more interesting. There the classical current re- 
mains indistinguishable from zero at all times, but the 
quantum acceleration is substantial. Note also that these 
"7;-classical" results have nothing to do with the true clas- 
sical limit of a DKRM, because here the DKRM is always 
on the main quantum resonance. Indeed, the ?7-classical 
Hamiltonian is given by Eq. ([^5]) . whereas the true clas- 
sical Hamiltonian of the DKRM should take exactly the 
same form as Eq. (1). 

We have also checked that if we choose = 0, tt in- 
stead, then both the classical and quantum acceleration 
seen in the examples in Fig. la do vanish. This confirms 
our previous discussion on a symmetry property of our 
new RA model. 



B. Dependence of Acceleration Rate on K and L 

To further explore the dynamical aspects of our new 
RA model, we have carried out detailed studies of how 
the ratchet acceleration rate depends on the system pa- 
rameters. The computational examples shown in Fig. 1 
motivate us to define the quantum current accelerate rate 
as follows: 

Rq EE d{p{t))/dt. (24) 

For the sake of comparison we also define the 77-classical 
current acceleration rate as 

i?e = d{p%t))/dt. (25) 

Computationally, these rates are determined as the aver- 
age linear increase rate over the time range 1000 < t/T < 
2000. Once the linear acceleration rates are obtained, we 
then check, in many cases, to see if the dynamics over a 
much longer time scale still accelerates the current with 
the same rate. Most often this is indeed the case, but 
some negative cases due to transient effects will be men- 
tioned below. 

Figure 2 shows the contour plots of \Rq\ and \Rc\ thus 
obtained as a function of K and L. A number of in- 
teresting features can be observed from Fig. 2. First, 
significant classical ratchet acceleration (see Fig. 2a) ex- 
ists only for those parameter regimes close to the K axis 
or the L axis. That is, at least one of the two values of K 
and L should be small for a considerable classical ratchet 
acceleration to emerge. But even that condition does not 
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FIG. 2: (Color online) (a): Absolute values of the classical 
current acceleration rate, denoted |-Rc|, as a function of K 
and L, for the RA model in the 77-classical limit, Eq. (|23|l . 
(b)-(d): Absolute values of the quantum current acceleration 
rate, denoted \Rq\, as a function of K and L, for the RA model 
in Eq. (I22j, with (b) ft = 1, (c) ?i = 27r/3, and (d) R = vr. In 
all the cases we set (j) = 7r/2. The contour scale lO"'^'^ is for 
values less than 10~^'^, the contour scale 10"^" is for values 
between lO"^'^ and 10"^ °, and so on. The contour scales 
used in panel (a) apply to other panels as well. 



FIG. 3: The phase space portrait of the »7-classical limit de- 
scribed by Eq. for 2L = A" = 1 in panel (a), 2J^ = L = 1 
in panel (b), 7^ = L = 1 in panel (c), and 21/ = isT = 4 in 
panel (d). For all the panels we choose (j> = 7r/2. In cases (a) 
and (b) phase space invariant curves extended in the momen- 
tum space can be clearly seen. In case (c) there is a web of 
separatrix structures and phase space invariant curves along 
the momentum space no longer exist. Case (d) represents 
a fully chaotic phase space if both K and L are sufficiently 
large. 



suffice. It is also clear from Fig. 2a that the regime of 
K ^ L should be excluded in order to have an appre- 
ciable \Rc\- The overall result is that in the parameter 
space defined by K and L, only a very small portion can 
yield considerable ratchet acceleration in the ?7-classical 
limit. Second, the quantum results shown in Fig. 2b- 
2d display many interest patterns. These patterns are 
absent in the classical case shown in Fig. 2a, and they 
vary strongly if we change h/n from an irrational value 
to a rational value. It can also be seen from Fig. 2b-2d 
that appreciable quantum ratchet acceleration occurs in 
a much larger parameter regime, often with \Rq\ > |-Rc|- 
Third, the quantum results share one feature with the 
classical result. That is, along the direction of K — L, 
\Rq\ is also seen to be small (typically much smaller than 
10-1-5). 

Some exceptions seem to be captured by Fig. 2b, where 
\Rq\ can become larger than lO^^'^ along the direction 
K = L. However, upon a careful investigation, we find 
that these exceptions are mainly caused by the particular 
way we numerically determine Rq. Indeed, if we follow 
the dynamics much longer (e.g. lO'* — 10® kicks), then the 
ratchet current tends to saturate for these exceptional 
cases, in contrast to the unbounded linear acceleration 
observed in other cases with K ^ L. Detailed investiga- 
tions of such transient effects in the ratchet acceleration 
are beyond the scope this work. The exact boundary be- 



tween bounded and unbounded quantum current acceler- 
ation can be an interesting and challenging mathematical 
problem. 

To shed more light on the results in Fig. 2, let us exam- 
ine in Fig. 3 the phase space structures of the ?7-classical 
limit of our RA model. The phase space structure of the 
entire phase space is just an infinite repetition (in both 
and p'^) of what is shown in Fig. 3. li K ^ L and either 
K or L is sufficiently small, then we always find phase 
space invariant curves extended in momentum. Trajec- 
tories moving along these invariant curves will display 
ballistic-like diffusion. Such a phase space feature dif- 
fers from that of the standard kicked Harper model. In 
the latter case the phase space invariant curves can lie 
parallel to the g'^-axis \i K <^ L. This difference is ex- 
pected, because the 77-classical Hamiltonian in Eq. pS)) 
is a twisted version of the kicked Harper model. 

Taking into account that q"^ = ±7r are equivalent points 
in the phase space, one can easily see that in both cases 
of Fig. 3a and Fig. 3b, there exist two bundles of phase 
space invariant curves, separated by a seperatrix struc- 
ture associated with some unstable fixed points. Remark- 
ably, the moving directions of the trajectories on the two 
bundles are opposite to each other. This feature is also 
consistent with the classical sum rule [2j . Based on these 
observations we are ready to explain the origin of the clas- 
sical accelerating ratchet current. In particular, because 
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the overlap of the p'^ = Q hne (the initial classical ensem- 
ble) with the two bundles of ballistic diffusion curves can 
be different, the effects of the two bundles of invariant 
curves cannot cancel out against each other and hence a 
net current develops. The current will increase linearly 
with time due to the ballistic nature of the phase space 
invariant curves extended in momentum. This under- 
standing is found to be consistent with an estimate using 
the intersection lengths between the p'^ = line and the 
two bundles of phase space invariant curves. This also 
clarifies the role of the parameter (p. As is evident from 
the expression of the 77-classical Hamiltonian in Eq. ((23)l . 
the net effect of the parameter is a shift of the phase 
space structure along the p'^-axis. So the parameter 
can be used to tune the unbalanced overlap of the initial 
ensemble with the two bundles of phase space invariant 
curves. 

By contrast, if = L (Fig. 3c), then before full chaos 
sets in, a web of separatrices emerge in the phase space 
and there are no longer phase space invariant curves ex- 
tended in momentum. As the value oi K = L increases, 
the chaotic layer associated with the web of separatri- 
ces becomes thicker and thicker. During this regular-to- 
chaos transition no phase space invariant curves extended 
in momentum will emerge. As a result, so long as K = L, 
ballistic-like diffusion cannot happen and a linear accel- 
eration of the ratchet current becomes impossible. This 
directly explains why the case oi K — L is so special for 
the T^-classical ratchet acceleration. Figure 2 also shows 
the absence of significant quantum ratchet acceleration 
in cases oi K = L. We believe that this quantum result 
is also due to the special classical phase space structure 
for k = L. 

Let us finally discuss the fully chaotic cases. One typ- 
ical example is shown in Fig. 3d. Analogous behavior 
can be found for other larger values of K and L. Clearly, 
phase space invariant curves are all broken in these fully 
chaotic cases. As a result, classical ratchet acceleration 
necessarily vanishes, as illustrated in Fig. lb. This ratio- 
nalizes the main message from Fig. 2a, i.e., appreciable 
classical ratchet acceleration exists only for a small frac- 
tion of the parameter space of K and L. 

Remarkably, in general full classical chaos does not for- 
bid quantum ratchet acceleration. As shown in Fig. 2b 
for a generic value of h (i.e., irrational with tt), large \Rq\ 
can be found for relatively large K and L, even when the 
associated ?7-classical dynamics becomes fully chaotic. A 
specific computational example has been shown in Fig. 
lb. This hence confirms a similar observation made in 
Ref. [l6j and constitutes another example of generic 
quantum violation of the classical sum rule 2j. In Ref. 
[16|, this violation was explained in terms of the con- 
centration of quantum amplitudes on the remanents of 
classical cantori-like phase space structures extended in 
momentum. Here, as K and L increase, the two bundles 
of phase space invariant curves will also generate classi- 
cal cantori-like structures. It can then expected that the 
classical cantori-like structures should be extended in the 




FIG. 4: (a) Dependence of the ratchet acceleration rate on 
the effective Planck constant h for = 7r/2 (dotted line) and 
(j) = 1 (solid line) . (b) Dependence of the ratchet acceleration 
rate on the system parameter tj), with h = n/2 (dash-dotted 
line) and h — 1 (solid line). In all the cases here 3L — K = 3h. 

momentum space for K < L as well as L < K. As such, 
the condition for quantum ratchet acceleration to occur 
in our new model should differ from that in Ref. p^ . 
Specifically, it is unnecessary to have a sufficiently large 
ratio K / L: a sufficiently small ratio K/L should also do. 
Given this understanding, the quantum ratchet acceler- 
ation here with full classical chaos is expected to display 
a "dual" symmetry in the K — L space. This is indeed 
the case: the pattern in Fig. 2b is symmetric along the 
line oi K — L. Note also that for non-generic values of 
Ti/tt, the quantum ratchet acceleration rates as shown in 
Fig. 2c and Fig. 2d can be even larger, despite the fully 
chaotic phase space in the underlying 77-classical limit. 



C. Dependence of Acceleration Rate on h and (fi 

Focusing on the quantum case, here we first examine 
the dependence of Rq on h. As already indicated by 
the drastic differences between Fig. 2b, Fig. 2c, and 
Fig. 2d, one might wonder if Rq is too sensitive to h 
such that experimental uncertainties in h may kill the 
ratchet acceleration altogether. To address this concern 
we show two typical computational results in Fig. 4a for 
two values of (j). In both cases some sharp peaks of Rq 
can be seen. These peaks are located at those values of h 
that are rational multiples of tt. Nevertheless, the overall 
?i-dependence of Rq does not show drastic oscillations. 
They can be varying smoothly over a considerably wide 
range of h. This feature also resembles what is found 
in Ref. [T6| . In the regime of very small fi, we have 
checked that Rq does approach the 77-classical accelera- 
tion rate i?c, thus establishing the expected quantum- 
classical correspondence. It should be stressed that the 
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?i-dependence of Rg shown here is a purely quantum ef- 
fect. To have a theory accounting for this ?i-dependence 
would be challenging but truly fascinating. 

A related question is whether or not the ratchet ac- 
celeration is robust to variations in the parameter cj) that 
characterizes the phase lag between two optical lattices of 
the same lattice constant. As demonstrated by the exam- 
ple shown in Fig. 4b, the 0-dependence of Rg is smooth 
in the entire range of cj). This is somewhat expected: in 
the ry-classical limit the parameter only shifts the phase 
space structure along the momentum axis. The conclu- 
sion is that small fluctuations in or ^ should not be a 
big concern for experimental studies. 

D. Effects of the Quasi-Momentum Spread in 
Cold-Atom Experiments 

In cold-atom experiments of kicked-rotor systems, the 
initial state cannot be exactly the state with zero quasi- 
momentum, even when the atoms are injected by a large- 
size Bose-Einstein condensate. Indeed, cold atoms are 
moving in real space, so their matter-wave state does not 
need to satisfy the periodic boundary condition inherent 
to a true kicked-rotor system. To motivate cold-atom ex- 
perimental studies of our new RA model, it becomes nec- 
essary to examine the detrimental effects of the nonzero 
quasi-momentum spread in the initial state. 

Because the kicking optical lattice potentials are al- 
ways periodic, the quasi-momentum of the cold atoms, 
denoted (3, is a conserved quantity [H, H^]. The dy- 
namics emanated from an initial state with a spread in 
(3 can then be easily simulated by considering each /3- 
component separately. To shed more light on this issue 
let us return to the DKRM propagator C/dkrm in Eq. ([2]). 
For each /? component, one now has p\n) = {n + /3)|n). 

This leads to the consequence that e~^'^^ ^ 1 under the 
quantum resonance condition Th = At:. Nevertheless, for 
Th = Att and for the potentials Vk (<?) and Vl (q) used in 
our RA model, it is enlightening to rewrite Eq. ^ as 

j-r/3 -i27r|^ ist -ikcoa{q+d,) -isi _i * cos{g) /f,c^ 

Except for the first factor, this expression is completely 
parallel to C/ra in Eq- (EH). Because the first factor 
is no longer unity and changes with /3, the first factor 
introduces de-phasing when a distribution of /3 values are 
averaged over. Therefore, it can be expected that the 
f3 spread will tend to saturate the accelerating ratchet 
current. 

Below we assume a Gaussian distribution of /3, with 
the variance denoted by A(3 and the mean value denoted 
by p. Results for a typical case are shown in Fig. 5. It 
is seen from Fig. 5a that a nonzero variance A/3 indeed 
induces the saturation of the quantum ratchet current. 
The exact saturation time increases, but slowly, with de- 
creasing A/3. For the parameters adopted in Fig. 5a, 
the typical saturation time is around 20 kicking periods 
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FIG. 5: Effects of the quasi-momentum spread on quantum 
ratchet currents with h = 1, 2L — K = 2, <j) = n/2. A{p{t)) = 
{p{t)) — (j5(0)). In panel (a), /3 = and the variance in /3 is 
given by A/3 = 0,0.0005,0.001,0.002 and 0.004 (from above 
to bottom). In panel (b), A(p(f)> at t = 7 with A/3 = 0.002 
(dashed line) and 0.01 (solid line) are shown as a function of 
the mean quasi-momentum /3 of the initial state. 

for A/3 ~ 0.002 (scaled by ?;.), a characteristic value of 
the P spread reported in a recent experiment using Bose- 
Einstein condensates [13] ■ Figure 5b also shows an inter- 
esting dependence of the ratchet current dX t — IT upon 
the mean quasi-momentum /3 of the initial state. The re- 
sult for A/3 = 0.01 is seen to be almost the same as that 
for A/3 = 0.002. This /3-dependence of the ratchet cur- 
rent at early times might be of interest to experimental 
studies as well. 

In future experiments the quasi-momentum spread can 
be made smaller than A/3 ~ 0.002 [23]. However, the re- 
sult in Fig. 5a indicates that even for A/3 ~ 0.0005, the 
saturation still sets in within a relatively short time scale. 
Interestingly, a previous theoretical RA model using a 
bichromatic optical lattice displays saturation at a sim- 
ilar time scale [9|]. This indicates that when faced with 
the detrimental effects of the quasi-momentum spread, 
the robustness of our new RA model without a bichro- 
matic optical lattice is similar to previous models with 
bichromatic optical lattices. 



IV. DOUBLE-KICKED ROTOR SYSTEMS ON 
HIGH-ORDER QUANTUM RESONANCES 

So far we have studied the ratchet transport in the 
DKRM under the main quantum resonance condition 
Th = Att. To motivate both theoretical and experimental 
studies in the future, in this short section we briefly dis- 
cuss an interesting extension of the current study. The 
extension is about the dynamics of ratchet current accel- 
eration in DKRM on high-order quantum resonances. In 
such extended cases, Th = At:v/ 11, where v and /i are two 
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FIG. 6: (Color online) Absolute values of the ratchet current 
acceleration rate, \Rq\ = \d{p)/dt\, for a double kicked rotor 
system at quantum anti-resonance TTi = 27r. The quantum 
propagator is given by Eq. (2). For the sake of comparison, 
the same parameter rescaling as in the cases of DKRM under 
the main quantum resonance is adopted. The rescaled param- 
eters are given by ?!, = 1 in panel (a), and ft = 27r/3 in panel 
(b). In both cases (j) = 7r/2. The meaning of the contour scale 
is the same as in Fig. 2, and the contour scales used in panel 
(a) apply to panel (b) as well. 

incommensurate integers. Under this high-order quan- 
tum resonance condition we find that analogous ratchet 
acceleration can be obtained as well, without using a 
bichromatic optical lattice. Thus, ratchet current ac- 
celeration itself may provide a useful tool for studies of 
high-order quantum resonances in DKRM. 

As an example in Fig. 6 we show \Rq \ under the quan- 
tum anti-resonance condition TTi = 27r, as a function of 
K and L. The kicking potentials VR-(g) and VL{q) are 
the same as those considered in Figs. 1-5. The case in 
Fig. 6a represents cases with a generic value of Ti, yield- 
ing appreciable \Rq\ in some regimes. In this case the 
detailed dependence of \Rq\ on K and L is seen to be 
rather complicated. The case of ?i — 27r/3 in Fig. 6b 
represents cases with nongeneric values of Ti. The associ- 
ated ratchet acceleration effect is seen to be larger over a 
wider regime. The dependence of li?^! ox\ K and L is also 
simpler than that seen in Fig. 6a. These results have no 
apparent connections with classical ratchet transport, in 
the sense that for high-order quantum resonances we can 
no longer define an 77-classical limit to guide our quali- 
tative understandings. However, as is evident from Fig. 
6, even in these quantum anti-resonance DKRM cases, 
a significant ratchet acceleration rate also requires the 
condition K ^ L. More studies of this extension will be 
carried out in the near future. 



V. CONCLUDING REMARKS 

To conclude, we have proposed and studied a new 
quantum ratchet accelerator model based on atom op- 



tics realizations of kicked-rotor systems. Unlike all pre- 
vious ratchet accelerator models, here we do not need to 
use a bichromatic optical lattice potential. Based on this 
advantage and the detailed computational studies pre- 
sented here, we believe that the cold-atom realization of 
our ratchet accelerator model is within the reach of to- 
day's state-of-the-art experiments [H, Indeed, the 
avenue of using atom optics to experimentally study a 
whole class of kicked-Harper-like models has been just 
opened up 0|, and the ratchet accelerator model pro- 
posed here seems to be a wonderful starting point along 
this direction. 

To have a linear acceleration of ratchet current in the 
77-classical limit of our model, we have shown that the 
phase space invariant curves extended in the momentum 
space are a necessary condition. Therefore, for kicking 
optical lattice potentials K cos{q) and Lcos{q + </>), an 
on-resonance double-kicked rotor with equal kicking am- 
plitudes K and L cannot yield an unbounded and lin- 
early increasing classical current. Instead, we need un- 
equal kicking amplitudes for accelerating and unbounded 
classical current to occur. This interesting requirement 
is also observed, but not fully explained, in the quan- 
tum dynamics. Given unequal kicking amplitudes, the 
quantum ratchet acceleration in our model can however 
persist for large kicking amplitudes, even when the rj- 
classical limit no longer has phase space invariant curves 
extended in the momentum space. This purely quantum 
effect is believed to be another example that remanents 
of classical phase space structures can dramatically im- 
pact the quantum dynamics. Considering these insights, 
we hope that our simple ratchet accelerator model will 
also motivate future theoretical work to better under- 
stand quantum transport and quantum-classical corre- 
spondence in classically chaotic systems. 
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